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Abstract 

The extremely large hierarchy observed in the fermion mass spectrum remains as one of the 
most puzzling and unresolved issues in particle physics. In a recent proposal, however, it was 
demonstrated that by introducing one Higgs doublet (or Private Higgs) per fermion this hierarchy 
could be made natural by making the Yukawa couplings between each fermion and its respective 
Higgs boson of order unity. Among the interesting predictions of the Private Higgs scenario is a 
variety of scalars which could be probed at future collider experiments and a possible dark matter 
candidate. In this paper, we study in some detail the dark matter sector of the Private Higgs 
model. We first calculate the annihilation cross sections of dark matter in this model and find 
that one can easily account for the observed density of dark matter in the Universe with relatively 
natural values of the model's parameters. Finally, we investigate the possibility of detecting Private 
Higgs dark matter indirectly via the observation of anomalous gamma rays originating from the 
galactic halo. We show that a substantial flux of photons can be produced from the annihilation 
of Private Higgs dark matter such that, if there is considerable clumping of dark matter in the 
galactic halo, the flux of these gamma rays could be observed by ground-based telescope arrays 
such as VERITAS and HESS. 
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I. INTRODUCTION 



One of the most puzzling issues in the Standard Model is the large hierarchy observed in 
the masses of fermions. For example, in the quark sector alone, the masses of the heaviest 
(top) and lightest (up) quarks are separated by nearly five orders of magnitude. Conversely, 
if one assumes that all fermions receive their mass via interactions with the same Higgs 
doublet (as in the Standard Model (SM)), the large hierarchy of masses observed in the 
fermion sector translates into a large hierarchy in the Yukawa couplings of the fermions. 

Recently, it has been proposed that the hierarchy of fermion masses can be made natural 
by extending the scalar sector of the SM to include one Higgs doublet (or Private Higgs 
(PH)) per fermion In this scenario, all of the Yukawa couplings can be made of 0(1) 
by tuning parameters of the model. In other words, the vacuum expectation values (vev's) 
of each respective PH field can be made to satisfy Vf ~ mf such that the hierarchy in the 
fermion mass spectrum becomes natural. 

The approach to electroweak symmetry breaking (EWSB) in the PH model is quite differ- 
ent than those of other multi-Higgs models. First, one introduces one gauge singlet scalar S q 
per quark flavor q and uses the vev's of these fields along with certain interactions between 
these fields and the various PH fields to induce "negative-mass-squared" instabilities. By 
using different terms in the Lagrangian for the top PH and non-top PH fields, one can easily 
explain the hierarchy in vev's by tuning certain parameters of the model. As a consequence 
of this approach, the lighter the fermion is, the heavier its associated PH particle must be 
in order to explain the smallness of the respective vev. In particular, the mass of the PH 
particle associated with the up quark can be shown to lie in the 10 2 — 10 3 TeV range which is 
definitely beyond the reach of the Large Hadron Collider (LHC). However, there is interest- 
ing phenomenology originating from the sector of the top and bottom PH fields along with 
the singlet scalars. In this work, we study a scenario where the physical spectrum of this 
sector contains a light SM-like Higgs boson, a heavy scalar Higgs boson, a pair of charged 
Higgs bosons and a pseudoscalar Higgs boson. The last three of these arise mainly from 
the bottom PH field and all have masses in the ~ TeV range. In addition to these, there 
are also two light scalars which are admixtures of the singlet states associated with the top 
and bottom quarks (S t and S&) . By construction, S t and are dark to interactions with 
SM gauge fields and fermions. While we will focus mainly on the light scalars in this work, 
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the heavier Higgs bosons could be probed at the LHC via production with bottom quarks 
(since the Yukawa coupling between the bottom quarks and the bottom PH field is of order 
unity) 1 . 

In order to avoid cross-talk between different quarks, the PH model contains a set of 
six discrete symmetries (one for each quark flavor). Under these symmetries, the right- 
handed quarks, their respective PH fields and the gauge singlet scalars S q are all odd, while 
all other SM fields are even. The existence of these discrete symmetries provides one of 
the most interesting features of the PH scenario which is the possibility of a dark matter 
(DM) candidate. Scalar DM was originally proposed over twenty years ago in Ref. [sj and 



las been studied more recently in several different scenarios including sing ! 
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et scalar DM 



151 ] . However, 



1 fl fl H, fl 3 and in the so-called Inert Doublet Model 0, Q, 
as we will demonstrate, the features of PHDM can be quite different from previously studied 
scenarios. 

The remainder of the paper is structured as follows. First, in Section [Til we review the 
structure of the PH model and demonstrate how E WSB is realized in this model. In section 

i— i n 

Section IIII[ utilizing the observations from WMAP [la ], we show that the PH model is 
able to account for all of the observed dark matter in the Universe for relatively natural 
values of the model's parameters. In addition, in Section HV] we consider the possibility 
of detecting PHDM via its annihilation into anomalous gamma rays in the galactic halo. 
We show that, with a favorable distribution of DM in the halo, PHDM could be detected 



by grounc 
telescope 



based telescopes, but is probably beyond the reach of the space-based GLAST 



13]. Finally, in Section [V] we conclude. 



II. THE MODEL 



The main goal of the Private Higgs model is to account for the extremely large hierarchy 
observed in the fermion mass spectrum [l|. For purposes of this paper, we will focus on 
the quark sector. In contrast to the SM, where one introduces a single scalar doublet 
which couples to all quarks, the PH scenario democratically introduces one Higgs doublet (f) q 
(q = u, d, s, c, t, b) per quark. All of the PH fields are assumed to have identical SU(2) xU(l) 

1 Presumably, the PH partner of the r could also provide interesting phenomenology; however, we will focus 
on the quark sector here. 
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quantum numbers as the SM Higgs. In addition to the PH fields, the scalar sector of the PH 
model also contains a set of gauge singlet scalars S q . In order to avoid cross talk between 
quarks of different flavors, a set of six discrete symmetries K q is imposed on the model. 
Under the K q symmetries, the right-handed quarks (U q , D q ) along with the PH fields and 
S q are all odd, i.e.: 

U q - —Uq (D q - -D q ) , <\> q - ~<\> q , S q - S q , (1) 

while all other fields are considered even. The Lagrangian which is symmetric under the K q 
symmetries is then given by: 

C = C S M-H-Y,i Y D H QL<t>DDq + Y^ H Q L 4> U U q ) 



-V{S q Aq), (2) 



where 4>u — i&2<f>u-> C-sm-h is the SM Lagrangian without the Higgs terms and Y£ H , Yjj H 
are Yukawa matrices. The scalar potential V(S q , <p q ) takes the form: 

v{s q A q ) = E{^^ + T^ + ^^ + ^ ( ^ )2 "^^ ff } 

+ a qq' S q S q' + 1 qq' S q S q' 'Pl'P q 1 + Xqq' Si 4>\<t> 



qq q q 1 iqq^q'^qYq'rq' i Aqq'^ q 'Y q Yq 

q+q' y 

+ aqq'^q'^q' + \q> < t>l < Pq4>\i < Pq> + Cqq'^q'^q'^q | + ^-C , 

where, for stability of the potential, a qq >,b qq r,c qq > < 0. In our analysis, we will assume these 
terms are small and neglect them in the following. 

In the PH model, instead of inducing EWSB through the usual "negative-mass-squared" 
approach where M| < 0, one utilizes the vev's of the singlet fields S q and the interactions 
between the S q s and the PH fields. In particular, for the top PH one assumes M| t > 
and induces EWSB through the g st and x q t couplings as well as the vev's of the S q fields. 
Thus, taking g st ,Xqt > and \M\ t - g st (S q ) 2 - J2q#Xqt{S q ) 2 = < 0, the top PH is 
forced to develop a negative-mass-squared instability which, in turn, spontaneously breaks 
the SU(2)l x U(l)y gauge symmetry. Therefore, in a sense, the top PH plays the role of 
the SM Higgs. 
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In general, the PH scenario can contain many new free parameters in addition to those of 
the SM. In order to simplify our analysis in the following sections, we will make a succession 
of approximations. Thus, our results will not probe the full parameter space of the PH 
scenario, but should be viewed as a first step in this direction. To begin, we follow Ref. [l| 
and assume that M? t <C Qst^ 2 which is in accordance with the symmetry breaking pattern 
discussed above. We also consider the case where g st ~ Xgt and , C 1. To give the S q 

potential in the S q — <p t sector reduces to: 



fields a vev, one introduces an instability Mj < such that, under our assumptions, the 



V(S q , <p t ) = ^ (s 2 - &*p\ 2 + X t (4<Pt? - 9stS 2 q 4<Pt , (3) 

where summation over quark flavor q is implicit and, in principle, the quantity vj is a bare 
parameter. Minimizing this potential, we find the conditions: 



dV(S q ,<j>- 



dS q 



X q s[( S ,) 2 -^ L )' 2 9st(4>t) 2 = 0, (4) 



q\2 



and: 

dV{S q ,<fH) 



Solving these equations for the individual vev's we find: 



= 2\ t (<p t ) 2 -g st (S q ) 2 = 0. (5) 

(S q ),{<t>t) 



2 \A|A t — g 2 st 
for the vev of S q and for the top PH vev: 



fa 9 )V 9stX q s \ vj 
4 \X 9 s K-9 2 stJ 2 



> 2 = ^U^F^r (7) 



Note that, for simplicity, we have identified the individual vev's vj with one common pa- 
rameter v s . Finally, we also note the relationship between v s and v^. 

2 9st 2 /o\ 

v h = ^v s • (8) 

Next, we consider the non-top PH fields which acquire their vev's in a slightly different 
manner. First, as in the case of the top PH, the mass parameter M| is assumed to be 
positive. However, for the (fi q fields (where q ^ t), one imposes the condition M| > g sq v 2 
in contrast to the case of the top PH. Then, vev's for the non-top PH fields are induced 
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through the cubic term 7^ and the vev's v s and Vh- Again, to simplify our analysis, we will 
make some assumptions. Specifically, we will assume that: 



Ml > g sq v 2 s , \ q (9) 



which is consistent with the symmetry breaking pattern discussed above. Then, after S q 
and (pt pick up vev's, the relevant part of the <p q potential is: 

Minimizing this potential, the vev's for the non-top Higgs fields are found to be: 

/ , \ Iqt V h V 2 s _ 

Eq. (TTTT) summarizes the main result of the PH scenario. By having the parameter 7^ to be 
small while keeping M| large, one is able to make all Yukawa couplings (which are given by 
m q /v q ) of 0(1) without fine-tuning. As a consequence of this relation, one can show from 
Eq. ( ITT]) that the lighter quarks have associated PH particles in the 10 2 — 10 3 TeV range 
which are definitely beyond the reach of current or future experiments pj. However, the 
masses from the <p t — S q sector can be naturally light (100's GeV), while the bottom PH 
particle can have masses in the TeV range. 

Finally, inserting Eqs. (J7J) and (TTTT) into the Lagrangian of Eq. ([2]), it is easy to show that 
the W mass is given in the PH model by: 

(12) 



Obviously, the leading term in the sum comes from the bottom PH; however, even in this 
case, the contribution is of order m 2 /m 2 ~ 0.001. Thus, the contributions to EWSB from 
quarks lighter than the top are negligible and our statement above that the role of the SM 
Higgs boson is being played by the top PH is verified. 



A. Mass Eigenstates and Their Interactions 

In this section, we study the top-bottom sector of the PH model in some detail. In 
particular, we will consider the case where A^. = X b s = A5 and Ag <^ A| for q ^ t,b. Thus, 
the gauge singlet scalars associated with the lighter quarks become heavy and effectively 
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decouple from our analysis. Then, under our assumptions, the scalar potential in the top- 
bottom sector reduces to: 



V(S q , <f> t , 4>b) 



4 



2\ 2 



s 2 t 



Si, 



,2\ 2- 



+ \ t (4^t) 2 + lMi4<j> b 



(13) 



- a-tb( s t s b) - ltbS t S b { 



fyt + 4>\<Pb) - 9st 



To begin, we expand the PH Higgs fields in the usual way: 



-L(v h + h t + i X °) 



(14) 




(15) 



while the singlet fields are expanded as: 



S, 



1 



t.b 



V2 



[V s + <?t,b) ■ 



(16) 



In the above expansions, ux^ and x° are assumed to play the roles of the usual Goldstone 
bosons which are eaten by the W and Z, while H and Af, are charged and pseudoscalar 
Higgs bosons, respectively. Both the H ± and Af, will have masses on the order of ~ 
TeV and could provide interesting phenomenology at the LHC Note that we are neglecting 
mixing between the "pure Goldstones" (cu 1 * 1 , x°) and the "physical Higgs bosons" (H ± ,Ab). 

2 

These mixings are typically of order and, thus, are extremely small. 

Inserting the expansions of the Higgs fields from Eqs. ( 1141 - ( |T6l) into ( |T3l) . we first extract 
the mass terms of the Goldstone bosons which we require to vanish: 

1 



x tv h - -g s tvt 



0. 



(17) 



Note that this equation is in agreement with Eq. ([8]). 



Next, we could attempt to diagonalize the full 4x4 mass matrix in the (ht, Hb,<Jt,<Jb) 
basis. However, as shown in Ref. |l|, f° r values of the model parameters that we consider 
in our analysis most of the mixings between the various scalars are negligible. In particular, 
the mixing between h t and a q is negligible provided: 



8g 3 st (2g st - X s ) 2 X t (for q = t) 



(18) 



and 



76t( — ) (for q = b). (19) 



Similarly, the mixing between a q and H b can be neglected provided: 

(20) 

All of these conditions are satisfied for the parameter choices in our analysis, hence we choose 
to neglect the above mixings. However, it should be noted that the mixing between ht and 
Hb serves to reproduce the small SM coupling between bottom quarks and the SM-like Higgs 
boson. In the following, we identify h t and H b with approximate mass eigenstates h° and 
H° respectively and assume the coupling between h° and a pair of bottom quarks takes its 
SM value. 

Finally, there can be substantial mixing between the two singlet scalars at and a b via 
the a s tb and jtb terms. Diagonalizing the 2x2 mass matrix in this sector, we find two mass 
eigenstates (Si and S 2 ) with mass eigenvalues: 



\ {>?s ~ 9stv 2 h - a s tb v^j - (^a s tb v 2 + -j= v hV^j sin 2a , (21) 
\ (4 - 9stv 2 h - a s tb v^ + (a s tb v 2 s + ^v h v b \ sin 2a . (22) 



Note that for a s tb ^tb > and < a < 7r/2, Si plays the role of the lightest PH parti- 
cle (LPHP) which is stable against decay and, thus, provides a candidate for DM. Using 
Eqs. (1211) and (122]) . we can exchange two of the free parameters (e.g., A5 and a s tb ) for the 
masses of the two singlet scalars. This is the approach we will take. Therefore, in the 
analysis to follow, we will take as our free parameters the masses ms n ms 2 as well as the 
couplings g st and 7^ and the mixing angle a. Note that the conditions for small mixings 
between the <T g 's and the PH fields forces g st to take small values. 

III. PRIVATE HIGGS DARK MATTER 

As mentioned earlier, one of the most interesting aspects of the PH scenario is the prospect 
of a Weakly Interacting Massive Particle (WIMP) with masses in the expected natural range 
for DM. In this context, the PH model is similar to other scalar DM models such as the 
gauge singlet models of Refs. Q, S Q, B, S fl 3 and the Inert Doublet Model (IDM) 




FIG. 1: Leading s-channel processes which maintain the singlet scalar Si in equilibrium with the 
rest of the cosmic fluid. 
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Q,B,Q. 



12, hj, hj, hj. In this section, we calculate the annihilation cross sections of PHDM 



into SM particles and show that, for relatively natural values of the model parameters, one 
can account for all of the observed dark matter in the Universe. In the next section, we 
investigate the possibility of indirectly detecting PHDM via its annihilation into anomalous 
gamma rays in the galactic halo. 

First, let us consider the present relic abundance of PHDM in the Universe. In the fol- 
lowing, we will assume that the mass splitting m^ 2 -m^ is large enough that coannihilation 
reactions between Si and £2 do not significantly affect the relic abundance. These effects 
will be considered in future work. In the early Universe, the singlet scalar Si would have 
been in equilibrium with the rest of the cosmic fluid. This equilibrium is maintained via £1 
pair-annihilation and pair-creation reactions which proceed through the s-channel exchange 
of the SM-like Higgs h°. The leading 2^2 s-channel reactions which contribute to these 
processes are shown in Fig. [TJ 

The present relic abundance of PHDM is determined by the pair-annihilation rates in the 
non-relativistic limit. The rates for each allowed channel are given in the non-relativistic 
limit as: 

a(X) = limo-(£i£i -> X)u (23) 

where u is the relative velocity of the annihilating particles. The total annihilation cross 
section is then given by summing over each of the allowed channels. Computing the cross 
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sections for the diagrams in Fig. [TJ we find: 



a{W + W ) = o ^; 2 !A _ 2 _ 2 A2 ; ^ 2 _ 2 ( 1 - n w + -n w ) , (24) 



2nvl 


(4m| a 








V 1 - V>z 




(4m| a 


- m l«) 2 + r2 h oKo 


■S'liSih 








(477^ 





.(zz> = ,,, ,, \\ J r2 , + (2r,) 



a(//) = 71 2 2 7 2 \_r2 2 . (26) 

where = m 2 / m| , T^o is the width of the h° for which we use SM values and the expression 
for the coupling <7ei£iA° is given by: 

gx^ho = -g st Vh + —f=-c a s a . (27) 



The WMAP collaboration 



provides a very precise determination of the present DM 



abundance which, at the two-sigma level, is given by: 

fl DM h 2 = 0.111 ±0.018. (28) 



As shown in Ref. [18j, for a generic model of DM, the present abundance of DM is mainly 
determined by Jo (the angular momentum of the dominant partial wave contributing to DM 
annihilation) and the annihilation cross section. In contrast, the relic abundance depends 
only weakly on the mass or spin of the DM particle. Thus, the very precise constraints from 
WMAP on Q^Mh 2 translate into very precise constraints on the quantity a = J2 x a(X) 
depending on the value of Jo- 111 particular, for an "s-wave annihilator" (J = 0) such as 
the case considered here, the WMAP measurement translates into the bounds: 

a = 0.8±0.1pb, (29) 



nearly independent of the mass or spin of the DM particle (see Fig. 1 of Ref. [181]). 

In Fig. [21 we plot the values of a(X) for two different values of the coupling g st and 
several different values of the mixing angle a. In these plots, we have set the bottom PH 
vev equal to the bottom quark mass and jtb = 1. The horizontal dashed lines indicate the 
limits on a(X) from Eq. ( J291) . Clearly, from Eqs. (I241 - (1271) . we see that the annihilation 
cross sections depend quadratically on g st . This is evident in the plot of Fig. [2] where we see 
a small shift in g st results in a large shift in the ranges of m^j allowed by the WMAP data. 
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Finally, we note that these plots are only meant to show that for relatively natural values 
of the model parameters it is indeed possible to account for the observed density of DM in 
the Universe. A full scan of the PH parameter space would probably find other choices of 
parameters which could fulfill the constraints from Eq. f[2§l) . 



IV. INDIRECT DETECTION OF PHDM 



Next, we would like to investigate the possibility of detecting PHDM. We will focus here 
on indirect detection and save an analysis of direct detection for future work. 

As we have seen, the annihilation rates for PHDM are approximately velocity-independent 
in the non-relativistic regime. In general, this implies that DM collected in galactic halos 
has a substantial probability to pair-annihilate resulting in anomalous high-energy cosmic 
rays which can be distinguished from astrophysical backgrounds. In particular, gamma rays 
from these annihilations provide a chance to extract information about DM, since they can 
travel over galactic scales without scattering. 

The production of gamma rays from E1S1 annihilation can originate from several different 
processes (including hadronization, factorization and radiation from final-state particles). 
However, for simplicity, we will assume the dominate source is from direct annihilation 
into a two-body final state as shown in Fig. [Hf. Note that, under our assumptions, only 
SM particles circulate the loop. In the full parameter space of the PH model, it would 
be possible to have charged Higgs circulating the loop. However, their couplings to h° are 
always of order v^/vh or v^/vg and, thus, can be safely ignored in comparison to the SM 
loops. The cross section for photon pair-production in the PH scenario can be written as: 

_ 2^ iSife0 f (ft -» 77 ) 

where the expression for g^^tfi is given above and, in the non-relativistic regime, s ~ 4m| . 
The hat on T indicates that one should replace m^o — > y/s in the standard expressions for 
on-shell Higgs decays. The expressions needed to construct T(h° — > 77) can be found in 



several reviews (e.g., see Ref. 19[]) and, hence, we will not repeat them here. 



2 Here, we concentrate on the dominant 77 signal and save a discussion of the Zj and/or h 0r y channels for 
future work. 
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FIG. 2: The annihilation cross section as a function of the Si mass and the mixing parameter 
a. The dashed horizontal lines indicate the WMAP constraints on the annihilation cross sections 
given by Eq. [29J. 



12 



FIG. 3: Diagrams which dominate photon pair-production in the Si annihilation in the galactic 
halo. 



Next, we would like to compute the flux of photons observed on Earth from Si annihi- 
lation in the galactic halo. The monochromatic flux due to the 77 final state, observed by 
a telescope with a line of sight parameterized by $ = (6, <fi) and a field of view AQ can be 
written as 



20(: 

$ = (1.1 x lO- 9 s-W 2 )(^) ^ 10 ° m GeV y^,AQ)AQ, (31) 

where the dependence of the flux on the halo dark matter density distribution is contained in 
J. Many models predict a large spike in the DM density in the neighborhood of the galactic 
center, making the line of sight towards the center of the galaxy the preferred one. However, 
the features of the peak are highly model-dependent resulting in values of J ranging from 
or AQ = 10~ 3 sr (typical for ground-based atmospheric Cerenkov telescopes) 



3 to 



21 



22 



7 



23 



24]. 



The monochromatic photon flux predicted from PHDM annihilation for the two values 
of g st studied previously are shown in Fig. HI For these plots, we have assumed there is no 
substantial spiking in the galactic center (i.e., JAQ = 1). In the energy range considered 
in these plots, ground-based atmospheric Cherenkov telescopes (such as VERITAS [25| and 



HESS 2g|) typically have a flux sensitivity down to the 10 -12 s _1 cm -2 level. On the other 
hand, the upcoming space-based telescope GLAST [l^] is limited by statistics to the 1CT 10 
s _1 cm -2 level over the energy range considered. From these plots, it is clear that without a 
substantial spike in the galactic center, PHDM will be difficult to observe in either ground- 
or space-based observatories. However, if the halo does exhibit a substantial spike or strong 
clumping (e.g., if J > 10 5 at AQ ~ 10~3), PHDM could be observed at ground-based 
telescopes assuming small values of g st and relatively light masses (ms a — 100 — 120 GeV). 
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FIG. 4: The flux of monochromatic photons from the reaction EiEi — * 77 for JAQ = 1 for two 
different values of the coupling g st . 
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V. CONCLUSIONS 



The Private Higgs model attempts to address the large hierarchy observed in the fermion 
mass spectrum by introducing one Higgs doublet for each fermion. EWSB is achieved not 
by the usual "negative-mass-squared" approach, but by introducing a set of gauge singlet 
scalars and using the vev's of these fields and their interactions with the PH fields to induce 
instabilities. In order to avoid cross-talk between quarks of different flavors, one also intro- 
duces a set of discrete symmetries. This provides one of the most interesting features of the 
Private Higgs model: a possible dark matter candidate. 

In this paper, we have begun an investigation of the PH dark matter sector. We found that 
for relatively natural values of the model's parameters the PH model provides a candidate 
which can account for the relic density of dark matter observed in the present Universe. 
To show this, we calculated the annihilation cross section for PHDM into SM particles and 
compared to limits on the cross section which can be obtained from the WMAP observations. 

Finally, we investigated the possibility of detecting PHDM via anomalous gamma rays 
originating from the annihilation of PHDM in the galactic halo. While the observation of 
these gamma rays may be difficult for the space-based GLAST observatory, we showed that 
evidence of PHDM could be observed at ground-based atmospheric Cerenkov telescopes such 
as VERITAS and HESS if there is substantial clustering of dark matter in the galactic halo. 
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